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0\ Abstract 
> . 

, The evaluation of hydrodynamic transport coefficients in relativistic field 

! theory, and the emergence of an effective kinetic theory description, is exam- 

^ ! 

ly-^ , ined. Even in a weakly-coupled scalar field theory, interesting subtleties arise 

at high temperatures where thermal renormalization effects are important. 

^ ', In this domain, a kinetic theory description in terms of the fundamental 

(D ! 

r~| , particles ceases to be valid, but one may derive an effective kinetic theory 



X 



describing excitations with temperature dependent properties. While the 



^ I shear viscosity depends on the elastic scattering of typical excitations whose 

kinetic energies are comparable to the temperature, the bulk viscosity is 
sensitive to particle non-conserving processes at small energies. As a re- 
sult, the shear and the bulk viscosities have very different dependence on 
the interaction strength and temperature, with the bulk viscosity provid- 
ing an especially sensitive test of the validity of an effective kinetic theory 
description. 
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I. INTRODUCTION 



In a weakly coupled quantum field theory, one would expect to be able to compute most 
physical observables starting from first principles. However, at sufficiently high tempera- 
tures, in even the simplest scalar field theory, the correct evaluation of transport coefficients 
characterizing long wavelength hydrodynamic behavior is quite subtle. Only recently has 
a thorough diagrammatic analysis of the bulk and shear viscosity appeared which is 
valid at temperatures where thermal renormalization effects are important. The purpose 
of this paper is to discuss the physical interpretation of the results of and to describe 
the formulation of an effective kinetic theory which properly incorporates thermal renor- 
malization effects and which generates the correct weak coupling behavior of both the bulk 
and shear viscosities. 

Existing literature in this area is somewhat sparse, particularly on aspects which are 
unique to relativistic quantum field theories. Consequently, we have tried to make the 
presentation reasonably self contained, and briefly review necessary background material. 
For simplicity, nearly all discussion will be limited to the case of a real scalar theory with 
cubic and quartic self-interactions, 

-C = \{d<pf + \mU' + l^^' + ^<P\ (1.1) 

with A <^ 1, mQ positive, and g"^ = 0{\'ml). Since the theory is weakly coupled, the physical 
(zero temperature) mass of the resulting scalar particles equals mo (after renormalization) 
up to radiative corrections, mphys = mo (1 + 0(A)). 

At non-zero temperature, the equilibrium state of this theory may be regarded as a 
fluid (or gas) of interacting spinless bosons. For fixed values of the coupling constants, 
the pressure, energy density, and other thermodynamic observables depend only on the 
temperature.^ It will be helpful to distinguish various ranges of temperature: 



Since the scalar field is real, particles are their own antiparticles and there is no conserved number 
operator or charge to which one could couple a chemical potential. 
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i) < T <^ fnpiiys. The system is non-relativistic and dilute. The equihbrium particle 
density is exponentially small, n ~ (mphysT)^/^e~™phy=/^. 

ii) mphys ^ T mphys/VX. The system is relativistic, but thermal corrections to the 
effective particle mass (or scattering amplitudes) are negligible. 

in) T ^ Tnpi^ys/ VX. The thermal correction to the particle mass, of order VXT, is 
comparable to the zero temperature mass. The system may no longer be regarded 
as a weakly interacting collection of the underlying fundamental particles. 

iv) mphys/ <^T <^ m-phys/ A. The zero temperature mass is negligible compared to the 
thermal mass shift, but the zero temperature mass still dominates the trace anomaly 
of/5(A)074! in (T/). 

v) T mphys/ A. The zero temperature mass is neghgible even in {T^^). 

The most interesting domains, from a theoretical perspective, are the high temperature 
ranges {iii-v) where thermal renormalization effects are important. Table | summarizes 
the qualitative behavior of various quantities at these temperatures. 



particle density 


n 


= 0(T3) 


energy density 


£ 


= 0(T4) 


effective particle mass 




= 0(Ai/2T) 


on-shell self energy 


S(p) 


= 0{XT^)+tO{X^T^) 


thermal width (p = O(mth)) 




= 0(A3/2T) 


thermal width {p = 0{T)) 




= 0{X^T) 


mean free time {p = 0{T)) 


Tf 


= 0(A-2T-i) 


elastic cross section {p = 0{T)) 


cr 


= 0(A2T-2) 


speed of sound 


^s 


= 1/VS + OiX) 


shear viscosity 


V 


= 0(A-2t3) 


bulk viscosity (T > 0(mphys/A)) 


C 


= 0{XT^ In^A) 


bulk viscosity (T = (mphys /\/A)) 


C 


= 0{ml^^,T-^X-^/^\n'X) 



TABLE I. Scaling behavior of various quantities in high temperature scalar field theory. The estimates 
hold in the domain T > ^riphys/VA where the one-loop thermal contribution dominates the (real part of 
the) single particle self energy S(p). If the scalar field has only quartic interactions, then the last result 
for the bulk viscosity acquires an additional factor of X^-^^'^. See section Mfor more detailed expressions. 
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The scaling behavior of the effective particle mass and thermal width will be essential 
ingredients in the following discussion. The thermal width is the inverse of the mean 
free time between scattering (up to statistical factors) and equals the displacement of the 
single particle pole away from the real frequency axis. The size of the thermal width follows 
directly from the imaginary part of the on-shell self energy divided by the particle energy. 
For weak coupling, the thermal width is small compared to the effective mass because the 
imaginary part of the on-shell self energy first arises from two-loop graphs, whereas the real 
part has one-loop contributions. The results displayed for the shear and bulk viscosities 
will be discussed in detail in section 0. 

II. TRANSPORT COEFFICIENTS AND BASIC KINETIC THEORY 

In a fluid with no conserved particle number, the stress-energy tensor T^^ is the only 
locally conserved current, and fluctuations in the energy and momentum densities are the 
only hydrodynamic modes.[| Two transport coefficients, the shear and bulk viscosities, 
(denoted 1] and (, respectively) characterize the resulting hydrodynamic response.^ If the 
system is slightly perturbed from equilibrium, then the non-equilibrium expectation of T^j^ 
will satisfy the constitutive relation (in a local fluid rest frame), 

= S^A-P) - (V.(n-) + V,(T°,) - |5.,V'(T°,)) - ^^%V'(T°,) , (2.1) 

together with the exact conservation law, d^T^^) = 0. Here Tij is the spatial part of 
the stress-energy tensor, e = Too is the energy density, and {V) is the local equilibrium 
pressure. The constitutive relation ( p.l|) is valid for small fluctuations in the limit in which 



^ These are fluctuations whose relaxation time diverges as the wavelength of the fluctuation increases. 
Such fluctuations determine the behavior of the system at arbitrarily long times and large distances. 

^Because there is no conserved particle number, thermal conductivity is not an independent transport 
coefficient. 
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the scale of the variation in {T^J) is arbitrarily large compared to microscopic length scales 
(such as the mean free path of excitations). 

The shear viscosity 7] characterizes the diffusive relaxation of transverse momentum 
density fluctuations; 7]/ (e+V) is the diffusion constant for such shear fluctuations. The bulk 
viscosity ( characterizes the departure from equilibrium during a uniform expansion. If the 
divergence of the fluid flow is constant, then the pressure differs from the local equilibrium 
value by the bulk viscosity times the expansion rate, or ^ V*(T°)/ {e+V). Both shear and 
bulk viscosity contribute to the attenuation of sound waves; the decay rate of sound waves 
with wavenumber k is k^(|r7+C) /{e+V) 0. The bulk viscosity vanishes identically in a 
scale invariant theory ||^ . This follows from the vanishing trace of the stress-energy tensor, 
Tf^^ = 0, in any scale invariant theory, and reflects the fact that a uniform dilation of an 
equilibrium distribution function remains in equilibrium (at a modified temperature) if the 
dispersion relation is scale invariant. 

Transport coefficients are proportional to the mean free path of the scattering processes 
responsible for relaxation of the associated hydrodynamic modes. The shear viscosity is 
proportional to the two body elastic scattering mean free path. In (scale non-invariant) 
relativistic theories, the bulk viscosity is proportional to the mean free path for particle 
number changing processes. This may be understood by noting that after a uniform ex- 
pansion, a change in the total number of particles is required in order to re-equilibrate at 
a different temperature.0 Note that decreasing the interaction strength will increase the 
mean free paths and thus normally increase transport coefficients. Consequently, the weak 
coupling expansion of viscosities typically begins with negative powers of the coupling. 

Most textbook discussions of the evaluation of transport coefficients (such as 0) begin 
by assuming the validity of a kinetic theory description of the interacting fluid. One argues 



*Under uniform expansion, a non-relativistic gas of molecules relaxes by converting internal energy (vi- 
brational or rotational) into kinetic energy. In contrast, a relativistic gas of structureless particles relaxes 
by converting rest-mass energy into kinetic energy. 
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that the system may be characterized by a distribution function f{x,p) giving the phase 
space probabihty density of the fundamental particles comprising the fluid. Although 
written as if it depends on an arbitrary four momentum, the distribution function is only 



defined for on-shell particles, for which = Ep = yp'^+m^Y^y^. The time dependence of 
the distribution function is governed by a Boltzmann equation, 

p^J. Q 



f{x,p) = \ j^^^ dTu^spl^fi f2 (I+/3) - (I+/2) /s fp) , (2.2) 



where dri2^3p is the differential transition rate for particles of momenta ki and ^2 to 
scatter into momenta and p, 



dT 



1 



12^3p 



2Ep 



T{p,ks;k2,k,)\ (2n)%k,+k2-ks-p) , (2.3) 



and fi = f{x,ki), fp = f{x,p). The collision term (or the right hand side of ( |2.2|) ) van- 
ishes when the distribution function is an equilibrium Bose distribution with an (inverse) 
temperature f3 and flow velocity m'^, or 

fP{x,p) = n{\u^p,\), (2.4) 

with n{E) the usual Bose distribution function at inverse temperature /3, 

n(E) = (e^^-l)"' . (2.5) 

To extract transport coefficients, it is sufficient to consider perturbations away from 
equilibrium which are arbitrarily small and slowly varying. Writing the distribution func- 
tion as a local equilibrium piece plus a non-equilibrium correction, 

fi^^p) = f'mi^^p) {1 - X{x,p) [1 + /;?.)(x,p)]} , (2.6) 

one may linearize the Boltzmann equation and expand in powers of Vm or V/5. After 
using the conservation relation, d^T^'^ = 0, to express time derivatives in terms of spatial 
gradients,^ one finds that (in the fluid rest frame at a particular point x) 



'And imposing the Landau-Lifshitz condition T'^'^u^, = T^^u^ to make the decomposition ( |2.6| ) unique 
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x{x,p) = P{x) A{x,p) V-u(x) + P{x) B{x,p) p-V(u(x)-p) — \V-u{x 



(2.7) 



where the coefficient B multiplying the shear in the flow satisfies the linear inhomogeneous 
integral equation 



PiPj - Ip^^ij = / '^ri2^3p (l+ni) (1+^2) ra3(l+rap)~^ 

2 J123 



Bi,{p) + Bijiks) - B,,{k2) - Bi,{ki) 



(2. 



with Bij{p) = B{p){piPj — ^Sij), and all quantities evaluated at the point x. The coefficient 
A multiplying the divergence of the fiow satisfies an analogous integral equation, 



12 2 / 2 I 2 \ 

3P - (P +'^phys) 



^ dTi2^3pil+ni)il+n2)n3il+np)-^ 
2 Jus 



X 



A{p) + A{ks) - A{k2) - A{k^) 



(2.9) 



with Vs = {dV / deY^"^ the (local equilibrium) speed of sound, together with the constraint 

d^p 







Epn{Ep){l + n{Ep))A{p). 



(2.10) 



Finally, inserting the distribution function into the kinetic theory stress-energy tensor 



T'^'ix) 



P^p'^fi.x.p) 



(27r)3Ep 

and comparing with the constitutive relation (|2.1|) yields 



(2.11) 



13 f d^p 



15 J (27r)3Ej 



p2n(Ep)(l + n(Ep))i?(p) 



and 



C = /3 



(2.12a) 



(2.12b) 



Hence, in this kinetic theory treatment, a quantitative evaluation of viscosities requires 



solving the linear integral equations ( |2.8| ) and ( p.9| ) and then computing the final momen- 
tum integrals in ( |2 . 1 2| ) . These results ( |2.9| - |2.12| ) for the viscosities can only trusted within 
the domain of validity of the underlying Boltzmann equation ( |2.2| ). Basic assumptions 
underlying kinetic theory which must hold include the following. 
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a) The collision time is negligible compared to the the mean free time between collisions 
of the fundamental particles. 

b) Between collisions, particles may be regarded as propagating classically with definite 
momentum and energy. 

c) The on-shell energy and momentum of particles in between collisions satisfy the zero 
temperature free particle dispersion relation, Ep — ■sjp'^+m'^^y^. 

When the system is non-relativistic, T -C mphys, the mean free time is exponentially 
large (compared to the Compton time h/nipbysC^) and the above assumptions are well 
satisfied. In the relativistic domain, T > mphys, the situation is more involved. In this 
regime, the density of particles scales as T^, a typical two-body elastic scattering cross 
section is cr ~ X'^/T'^, and so the mean free time is ©(l/A^T) (or ©(l/A^/^T) for soft 
particles due to Bose-enhanced stimulated emission). In contrast, the typical collision time 
(determined by the variation of the phase shift with energy) is 0(A^/T); hence condition (a) 
is satisfied as long as the theory is weakly coupled. Quantum uncertainties in the energy or 
momentum of a particle propagating between collisions are negligible provided the kinetic 
energy times the mean free time is large compared to h. For particles with typical 0{T) 
energies, this condition again merely requires A <^ 1. But since the mean free time becomes 
arbitrarily small as the temperature increases, "soft" particles with momentum of order of 
their rest mass cannot be viewed as propagating classically when T > "niphys/A^. Moreover, 
standard kinetic theory fails long before this temperature is reached due to condition (c). 
The collision term in the Boltzmann equation summarizes the effects of scatterings in 
which particles change their momenta in a near-random manner which may be regarded 
as destroying phase coherence. It does not describe the coherent change in phase caused 
by exactly forward scattering. The amplitude for a soft particle to propagate through the 
surrounding medium will be modified due to phase shifts arising from forward scattering 
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interactions, and this will change the dispersion relation from the zero temperature form.|^ 
For a hot scalar theory, this is precisely the origin of the well-known thermal correction to 
the effective particle mass, 

-th(T)^ = mj,^, + ^ X (l +0(^)) . (2.13) 

For simplicity, contributions arising from the cubic coupling have not been displayed. Here 
(and henceforth) A and are renormalized couplings evaluated at the scale T.[] The 
thermal mass correction is negligible when T < O(mphys), but when T > 0(mphys/"\/A) 
the mass correction is significant and the standard Boltzmann equation fails to describe 
correctly the propagation of particles with soft O(mphys) momenta. It is important to note 
that forward scattering effects will also change the effective cross sections of soft particles 
propagating through the medium. 

When T ^ TTiphys, one might expect the inapplicability of the Boltzmann equation for 
soft particles to be irrelevant, since particles with O(mphys) momenta comprise a small 
0((mphys/T)^) fraction of the total. This is true for some physical observables (includ- 
ing thermodynamic quantities such as the pressure or energy density, and also the shear 
viscosity). However, as will be seen explicitly below, the bulk viscosity is predominately 
sensitive to soft O(mphys) momenta.^ 

There is an additional problem with the kinetic theory treatment for the bulk viscos- 
ity. The integral equation (|2.9|) has no solution! The kernel of the equation has a zero 
mode (which is not orthogonal to the source term). The zero mode is a consequence of 
the conservation of particle number. The original gip'^+Xip'^ theory does not have a con- 



^This, of course, is exactly how the index of refraction for hght is generated. 

renormahzation point of order T is needed to avoid large logarithms in higher order corrections. 
^This difference between the shear and bulk viscosity is easy to see from equations (p^- 2.1^ ). The factors 



of in the shear viscosity integrand ( 2.12a ) and the inhomogeneous term in ( |2.8| ) combine to suppress 
the contribution of soft momenta by four powers of |p| relative to the case of the bulk viscosity. 
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served particle number, but the Boltzmann equation (|2.2|) with only two-particle scattering 
terms included obviously does conserve the number of particles. As stated earlier, the bulk 
viscosity, which characterizes the relaxation of the system after a uniform expansion, is 
directly sensitive to particle number changing processes since a (scale non-invariant) sys- 
tem undergoing uniform expansion cannot re-equilibrate without changing the number of 
particles.0 Consequently, higher order particle number changing terms must be included 
in the Boltzmann equation ( p. 21) even though they are suppressed by additional powers of 
A. This will be described more explicitly below. 

In summary, standard kinetic theory (with number changing processes included) is 
adequate for calculating transport coefficients in a weakly coupled theory in the tem- 
perature regimes where T <^ "rn-p^ys/ ^/^i but not in the high temperature regimes with 
T > mphys/ V^- In order to derive transport coefficients in this domain, one should start 
directly from the underlying field theory. 

III. DIAGRAMMATIC EVALUATION OF TRANSPORT COEFFICIENTS 

The shear and bulk viscosities may be extracted from the zero momentum, small fre- 
quency limit of the spectral density of the equilibrium stress tensor-stress tensor correlation 
function. One finds that [IH 

( dSe'"*([7ri^(t,x),7ri^(0)])eq , (3.1a) 
d^a;e*-*([P(t,x),P(0)])eq. (3.1b) 

^If the number of particles is conserved (as in a non-relativistic field theory), equilibrium states will 
depend on a chemical potential as well as the temperature, and the zero mode in ( |2.9| ) will be removed by 
the additional subsidiary condition on x needed to make the local chemical potential uniquely defined. In 
such a theory, a uniform expansion will, of course, produce a change in the chemical potential. 



ri = — hm — 

' 20 ^--0 uj 



and 



C = - lim — 

2 "^^0 u 
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Here, nim = Vz0 Vm0 — If^^mlV^)^ is the traceless part of the stress tensor and V = V—v'^e 
is the pressure minus the energy density times the square of the speed of sound.0 These 
Kubo relations provide the natural starting point for a field theory evaluation. The spec- 
tral density equals the discontinuity in the (Fourier transformed) stress-stress correlation 
function and has a perturbative expansion generated by the sum of cut diagrams with two 
insertions of T^j, @]. Naively, one would expect the leading weak coupling contribution to 
arise solely from the single one- loop diagram shown in Fig. 1. 



FIG. 1. The cut one loop diagram contribution to the viscosity. 

This is correct for generic values of the external 4-momentum, but is completely incorrect 
in the limit of vanishing external momentum and frequency. Finite temperature propaga- 
tors have poles (in frequency) with both +ie and —ie prescriptions. When the external 
4-momentum vanishes, the product of propagators corresponding to the graph in Fig. 1 
contains terms in which the contour of the frequency integration is pinched between coa- 
lescing poles, thereby producing an on-shell divergence. As always, such divergences have 
a simple physical origin [^. When a small momentum is introduced by an insertion of T^iy, 
an on-shell (bare) particle in the thermal medium can absorb the external momentum and 
become slightly off-shell. The amplitude is proportional to the length of time the parti- 
cle can remain off-shell. As the external momentum vanishes, the virtual particle moves 



-'^'^The Kubo relation (3.1b) is equally correct if the pressure V is used in place of "P, since commutators of 
the energy density vanish at zero momentum (due to energy conservation). However, as shown in Q and 
explained below, the particular choice of V given is appropriate for deriving an "effective" kinetic theory 
description. 
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on-shell and the integral over the propagation time diverges. 

However, at non-zero temperature, no excitation can actually propagate indefinitely 
through the thermal medium without suffering collisions off other excitations. In a scalar 
field theory, a single particle excitation of momentum k acquires a finite lifetime r^, or 
non-zero thermal width = l/rjt, due to the O(A^) imaginary part of the on-shell two- 
loop self-energy. To examine the limit of vanishing external momentum, one must resum 
the single particle self-energy insertions which will shift the poles in the single particle 
propagator from ±E^±ie to ±E]^±iTk (where E]^ = y^k^+m(T)^). This serves to 
regulate the apparent on-shell singularity, and makes the one-loop diagram in Fig. 1 yield 
a finite result proportional to the single particle lifetime. However, since the lifetime is 
0(1/A^) (for particles with 0{T) momenta) this means that higher loop diagrams can be 
just as important as the one-loop contribution if they are sufficiently infrared sensitive. 




FIG. 2. A typical cut ladder diagram for the shear viscosity in gip^+Xcf)'^ theory containing O(A^), 
0{g^X), and 0{g'^) "rungs". 

For the shear viscosity, a careful analysis shows that one must sum all cut "ladder- 
like" diagrams of the type illustrated in Fig. 2. See [|1| for details. This is similar to the 
situation in non-relativistic systems except that instead of having ladders built from 
an instantaneous two body interaction, one must deal with ladder graphs containing far 
more complicated "rungs" . Nevertheless, one may formally sum all cut ladder-like graphs 
by introducing an effective vertex V^^lk, q—k) satisfying a linear equation of the form 

V^{k, q-k) = I^{k, q-k) + j -^M{k-p) J^{p, q-p)V^{p, q-p) . (3.2) 

The effective vertex actually has a four components (in order to represent the four different 
choices for which legs are above and below the cut), while A4{k—p) and J^{p,q—p) are 
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4x4 matrices representing the rungs and side-rails of the ladder, respectively. These 
matrices have entries consisting of various products of cut and uncut propagators. The 
inhomogeneous term X7r(/c, q—k) represents the vertex factors corresponding to an insertion 
of the traceless stress tensor. The explicit form of each of these quantities may be found 
in Closing the two legs of the effective vertex with a second insertion of the traceless 
stress tensor produces the sum of all ladder-like graphs contributing to the shear viscosity, 
so that 

7/ = ^ lim hm / T.(fc, q-kyj'ik, q-k) V^{k, q-k) x (1 + 0(v^)) . (3.3) 

In the limit of vanishing external momentum q, one may perform the frequency integra- 
tion and extract the leading order behavior from the nearly pinching-pole contributions [|l| . 
Moreover, by using the finite temperature optical theorem |]l],|^ the 4x4 kernel AiJ^ may 
be shown to equal a rank one matrix (up to corrections subleading in A), thereby allow- 
ing one to reduce the equation to a single component, three dimensional integral equation. 
The result is identical to equation (p.8|) for the spin-two part of the Boltzmann distribution 



function, and (|2.12a|) for the kinetic theory shear viscosity provided one: 



a) identifies the shear response B{k) with the effective vertex divided by the imaginary 
part of the single particle self-energy, 

b) uses the thermal mass mth(T) instead of the zero-temperature mass in the dispersion 
relation defining on- shell momenta, and 

c) uses an effective temperature-dependent "scattering amplitude" equal to the usual 
tree-level amplitude but evaluated with finite temperature retarded propagators,Q 

^(Pi,P2;P3,P4) = >^-f {.Gr{.Pi+P2) + Gr{.Pi-Pa) + Gr{.Pi-P^)) , (3.4) 



"'^'^ Since the intermediate propagator in (3^) cannot go on-shell, the ie prescription in the retarded prop- 
agator is actually irrelevant. Note however, that using the real time Feynman propagator in the scattering 
amplitude is incorrect as this differs (off-shell) by a Bose distribution factor. 
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where Gji{pi+p2) = [{pi+P2)'^ + ^th]~^ g = g + \{(f)) is the "shifted" cubic 
couphng constant that resuhs when one shifts the field by its thermal expectation 
value {(p) in order to remove tadpole diagrams. 

The calculation of bulk viscosity requires considerably more care than the shear vis- 
cosity. In addition to ladder diagrams of the type shown in Fig. 2, one must also sum 
diagrams containing iterations of higher order number- changing scattering processes, and 
include thermal "vertex renormalization" subgraphs [Q]. Examples are shown in Fig. 3. 



FIG. 3. A typical graph containing 0{g^X^) and 0{g^) two-to-three particle "rungs", plus "thermal 
renormalization" of the stress tensor vertices. Graphs such as this contribute to the leading order weak 
coupling behavior of the bulk viscosity. 

Nevertheless, one may again sum all the relevant diagrams by introducing an effective 
vertex T>p{k, q—k) satisfying a linear equation of the same form as in ( ^3.21 ). The appropri- 
ate kernel now contains the previous O(A^) subdiagrams plus 0{g'^X^) number changing 
subdiagrams.Q The inhomogeneous term receives 0(A) corrections involving the one- loop 
contributions to the thermal mass and the speed of sound. These vertex corrections cannot 
be neglected at high temperatures because the speed of sound (squared) approaches 1/3, 
producing a cancellation in the leading O(p^) part of the inhomogeneous term ( |2.9| ). Con- 
sequently, an insertion of V (or pressure minus fg times the energy density) is 0(mpjjyj, 



Scattering amplitudes, strictly speaking, do not exist at non-zero temperature, since all excitations have 
finite lifetimes. However, in this weakly coupled theory, the effective scattering amplitude (3.4) provides 
a meaningful characterization of scattering processes which occur on time scales short compared to the 
single particle lifetime. 
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Or, for a pure X<p'^ theory, 0{X'^) two-to-four particle subdiagrams. 
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even when the loop momentum is of 0{T). Hence, vertex corrections which are O(AT^) 
can be comparable to the zeroth order term. 

Once again, one may perform the frequency integrations and extract the leading be- 
havior from the nearly pinching-pole contributions,!^ show that the resulting kernel is 
dominated by a rank one matrix, and reduce the equation to a single component, three 
dimensional integral equation. This has the same form as Eq. ( p.9|) (with A{p) identified 
with the effective vertex divided by the imaginary part of the self-energy) except that: 

a) In addition to the two particle elastic scattering term, the right hand side now con- 
tains a particle number changing term proportional to the square of the tree level 
two-to-three particle "scattering amplitude" ,Q 

Taat = iXg J2 GniPi+Pj) - if GR{pi+Pj)Gnipi+Pm) , (3.5) 

where again g is the shifted cubic coupling constant. 

b) The thermal mass is used in the dispersion relation for on-shell momenta, and in 
the retarded propagators appearing in the "effective" thermal scattering amplitudes 
(0) and (pD. 

c) The physical mass (squared) appearing in the source term is replaced by 

~ 2 _ 2 rTi2 ^^Iyi /o r\ 

m =m^^-T ——. (3.6) 



^^Subleading non-pinching pole terms in the kernel can be neglecting only if the inhomogeneous term is 
orthogonal to the zero modes of the reduced pinching-pole kernel (as well as orthogonal to the zero modes 
of the full kernel). Imposing this condition forces the energy density coefficient in the source V = V—VgE 
to equal the speed of sound (including one-loop corrections) 0. 

^^Or, for a pure A^'' theory, the two-to-four particle amplitude Tan — GB.{Pi+Pj+Pk)- Here 

all 6 momenta involved in the two-to-four scattering are regarded as incoming, and the sum runs over 10 



distinct partitions of the six momenta into two groups of three momenta. Similarly, the sums in Eq. (3.5) 
run over partitions of the five momenta into sets of 2 and 3 momenta, or 2, 2 and 1 momenta, respectively. 
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The "subtracted" mass rh'^ is a measure of the departure from scale invariance. The 
subtraction cancels the leading temperature dependence in m^j^, so that differs negligibly 
from mpjjyg when T < mphys/V^, and approaches rripj^yg — ^g'^/X for mphys/v^ <^ T <^ 
"^phys/A. At asymptotically large temperatures, T ^ mphys/A, the running of the quartic 
coupling in ( |2.13| ) dominates and = /5(A)T^/48, up to 0{VX) corrections. 
The resulting equation for the spin-0 response is 



2 I ~2\ 



(p +m 



E„ 



2 J 123 
X 

E. 



dTi2^3p (1+ni) (1+722) ns (1+% 



A{p) + Aih) - Aih) - Aih 

E f 

4 J 1234 

+ A{k^) + A(A;3) - A{k2) - A{h) 



X 



+ ^ / C?ri23«4p (1+ni) (1+^2) (l+ng) n4 (l+np)-^ 
J 1234 



X 



A{p) + A{k,) - ^(^3) - ^(A:2) - A{h] 



(3.7) 



with 



12^34p 



2E, 



1 |2 ^ d^k,' 



{2^mp^:-Po:i) , (3.8) 



(27r)3(2E,,; 

eic. In the pure quartic theory, the 2^3 particle terms are replaced by the corresponding 
2 4 particle contributions. Closing the effective vertex with an insertion of V yields the 
bulk viscosity. 



C = P 



|p2 - v', ip'+m')) n{Ep) (1 + n{Ep)) A{p) 



(3.9) 



which differs from ( p.l2b| ) by the replacement of rn^Yiys by 



^^The solution of ( |3.7|) for A{p) is only unique up to the addition of a zero mode contribution proportional 
to Ep. This has no effect on the bulk viscosity ( ^j| ) because the speed of sound satisfies the identity 

= j (fp/{2^f (ip2 - vl (p2+m2)) n{Ep) {I + n{Ep)) . (3.10) 

Nevertheless, the ambiguity in Aija) may be eliminated by imposing the Landau-Lifshitz condition for the 
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IV. EFFECTIVE KINETIC THEORY 



Before discussing the solutions of these hnearized equations for the hydrodynamic re- 
sponse, we wish to show how one may construct an effective kinetic theory for quasi-particle 
excitations which reproduces, at arbitrary temperature in a weakly coupled theory, the cor- 
rect hydrodynamic response. As usual, the quasi-particle distribution function f{x,p) will 
depend on an on-shell four-momentum p, but now the quasi-particle energy = Ep will 
be a function of both the spatial momentum p and an effective mass m{q), which in turn 
depends on a spacetime-dependent auxiliary field q{x): 



The auxiliary field q characterizes the effect of the forward scattering of a quasi-particle off 
other excitations in the medium, and depends self-consistently on the distribution function. 



This is just a non-equilibrium generalization of the usual thermal contribution to the scalar 
field propagator at coincident points, The quasi-particle Boltzmann equation can 

be written as 



The dispersion relation ( [4.1|) implies that dEp/dp = p/Ep and dEp/d:s. = {m/ Ep)'Vm. 
Hence, the spatial gradient of the effective mass acts like an external force which changes 
the momentum of propagating excitations. The collision term on the right hand side is the 
usual Boltzmann collision term with both 2 2 and 2 ^ 3 (or 2 ^ 4 for a pure quartic 
theory) particle processes included. 



effective theory described below. This reduces to the constraint 






(4.2) 




(4.3) 




(3.11) 
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AT{x,p) = \ j^^^ rfri2^3p (/l /2 (I + /3) + - (I + /1) (I+/2) h fp) 

+ I I rfri2^34p (h f2 (I + /3) (I+/4) - (I + /1) (I + /2) fs U fp) 

J1234 V / 

+ |/^^/ri23^4p(/l/2/3(l + /4) (l+fp) - (I+/1) (I+/2) (l+/3)/4/p). (4.4) 

The transition rates (for a given spacetime location x) are given by the usual definitions 
( p.3|) and ( p.8|) , with effective scattering amplitudes ( |3.4|) and ( |3.5|) computed using retarded 
free propagators containing the effective mass m{q{x)). 

This effective Boltzmann equation is to be combined with a modified definition of the 
kinetic theory stress-energy tensor, 

T^^ix) = (I ^0-p^p'^f{x,p))-g^^''U{q{x)). (4.5) 

A short exercise shows that the modified stress-energy tensor (|4.5|) is conserved provided 
the interaction energy U{q) satisfies dU/dq = —\q{dw?/dq), or 

U{q) = \ d(i (m2(g') - m'iq)) . (4.6) 

This is also the necessary consistency condition for ensuring that the variation of the total 
energy density with respect to the quasi-particle density yields the correct quasi-particle 
energy, Ep{x) = 6T^^{x) j 6f{x,p) 0, and in equilibrium, that the pressure satisfy the 
correct thermodynamic identity T{dV/dT) = e + V. 

The final ingredient needed to complete the definition of the effective kinetic theory is 
the dependence of the effective mass on the auxiliary field q. This is completely determined 
by the dependence of the equilibrium thermal mass m^h on the one-loop "bubble" (0(a;)^). 
In the pure quartic scalar theory, the thermal mass has the simple form, 

m^{q) =ml + l\q (4.7) 

(up to corrections suppressed by powers of A), while if cubic interactions are present one 
must first self-consistently expand the field about its thermal expectation value c = (0), 
leading to 
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m'^iq) = ml + gc+ lX{c^+q) , (4.8) 

with = mgC + \g{c^+q) + | Ac(c^+3g). As always, the couphng constants appearing 
in ( |4.7| ) and ( |4.8|) should be evaluated at a scale appropriate to the physics under con- 
sideration; the running of the quartic coupling affects even leading order results when 
q > mpjjyg/A^. In equilibrium, q ~ when T ^ mphys- Hence, the appropriate gener- 

alization is to regard the coupling as an implicit function of q satisfying (when q ^ ^phyj 

q—^ll3(X) = lboX' + OiX'). (4.9) 

The resulting effective mass in, for example, the massless pure quartic theory is 

where A = the renormalization group invariant scale of massless 0^ theory. 

This effective kinetic theory provides a consistent description of the non-equilibrium 
dynamics of a weakly coupled scalar field theory, including the propagation of slowly moving 
excitations, even when the effective mass of the excitations differs substantially from the 
zero-temperature mass, or varies significantly in space or time. Expanding about a local 
equilibrium distribution, as in (PTB|), and evaluating the effective stress energy tensor ( [4. 51 ) 
(carefully keeping track of the implicit dependence on the distribution function hiding in 
every factor of energy), leads to the fairly simple result 

T'^^{x)=T!f^^{x)- J j0^^n{E,){l+n{E,))x{x,p)(p''p^-uVT'^y (4.11) 

where T^^ = u^u" (s+V) + g^^ V is the local-equilibrium contribution. Expressing xi^iP) 
in terms of the shear and bulk amplitudes (c./. Eq. ( p.7|) ), and linearizing the effective 
Boltzmann equation in the hydrodynamic limit, yields exactly the same equations obtained 
in the previous section for the amplitudes v4(x, p) and B{x, p). When inserted into the stress 
tensor ( [4.11| ) one precisely obtains the previous results ( p.l2a| ) and ( |3.9[ ) for the shear and 



bulk viscosities. 
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V. RESULTS FOR VISCOSITIES 



Computing the bulk viscosity requires solving the integral equation ( p.7|) . Unlike the 
case of the shear viscosity, solving this equation is trivial because the kernel has a single 
small eigenvalue which is only displaced from zero due the inclusion of number changing 
processes. Hence, the solution is dominated by the projection onto the near-zero mode, 
leading to 

A{p) = -^{l-aEp) , (5.1) 



r 



AN 



where 



JWEp 

with Ip{p) = |p^— fg (p^+m^) the same source term as in Eq. ( p.TD , and Faat the total 
3 — »• 2 particle (or 4—^2 for pure A</)^) thermal reaction rate per unit volume. 



/ ^(Ep) ^v{p) , (5.2) 



X ([l+n{E,)] [l+n{E2)] n{E,) n{E^) n{E,) ) . (5.3) 



The constant a in ( [5.1D is undetermined by (|3.7D, but may be adjusted to satisfy ( |3.11| ). 



The bulk viscosity obtained by inserting into ( p.9|) is simply 



F2 

C = Pt^- (5.4) 
i / 



The final evaluation of the shear viscosity requires a numerical solution of the integral 
equation (2^) and the final integral ( 2.12a|) , while the bulk viscosity requires performing 



the rather involved phase space integral (|5.3| ) for the particle number changing reaction 
rate. Details of this evaluation may be found in [Q. 

Despite the need to resum self-energy insertions in order to cut off singularities in the 
original diagrams, the introduction of thermal corrections in the dispersion relation and 
scattering amplitude is actually irrelevant for the leading behavior of the shear viscosity 
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because the integrals ( |2.^ ) and ( p.l2a| ) are dominated by momenta of order T. This, 
however, is not the case for the bulk viscosity. At high temperature, the number changing 
reaction rate scales as 0{g'^\^T^ /m^]^ for the g(j)^+\(j)^ theory, and 0(A^T^/m^jj) for the 
pure theory due to its infrared sensitivity to the region where all momenta are O(mth)- 
The factor F appearing in the numerator is a measure of the violation of scale invariance of 
the theory, and behaves as 0(m^T^ ln(T/mth)) when mphys <^ Hence, the shear and 
bulk viscosities have very different behaviors throughout the high temperature region.|3 In 
pure theory. 



2^3 

r7 = a— X 



0(yA) + 0(mphys/T) 



(5.5) 



while 



T 



(5.6) 



when mphys <^ mphys/ A, and 



C = cAln2(7A) x 



0(VA) + 0(mphys/AT) 



(5.7) 



when T ^ mphys/A. The forms ( |5.5| ) and ( p.TD remain valid if cubic interactions are 
present, but the bulk viscosity in the intermediate regime mphys ^ T ^ mphys /A acquires 
dependence on the relative strength of cubic and quartic couplings. 



"One finds vl ^ \ - ^ to^/tt^T^ and F = -{m^T^ /6Tr^) [ln(2T/mth) - ^C(3)/7r2], ^j^g^ evaluating 
(3.10) and (5.2) for T 3> mphys, up to corrections suppressed by (VX) or mphys/T. 



"'^^They are also very different at low temperature. When T ^ mphys, the shear viscosity behaves like 
rj ~ mpj^yg(T/mphys)^^^/A^, but the bulk viscosity diverges exponentially as C gSmphys/T m^^^^/ X'^T'^ for 
pure A(/)* theory, or e™p''y=/-^ (mphys/T)^/^ mpj^y^/A^g^T for g(j)^+\(t>'^ theory. This is the bulk viscosity 
characterizing asymptotically long wavelength hydrodynamic fluctuations, appropriate for distances large 
compared to the mean free path for particle number changing interactions (which displays the same 
exponential divergence). Ordinary non-relativistic hydrodynamics (with a conserved particle number) is 
valid at distances small compared to this number changing mean free path but large compared to the elastic 
mean free path. It is, of course, this region and not the strict asymptotic domain which has practical utility. 
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C = <^)P^ln^(^)x [l + 0(v^) + 0K,,3/T) + 0(AT/mp,,3)], (5.8) 

with d{x) a non-trivial dimensionless function. 

A numerical evaluation of equation (|2.8|) , (|2.12a| ), (|5^ ), and for the pure quartic 
theory yields the values 

a = 3.04 X 10^ , (5.9a) 

6 = 5.5xl0^ (5.9b) 

c = 6/(6 (327r)^) = 8.9 X 10-^ (5.9c) 

gi5C(3)/2.2 ^2 = 1.2465, (5.9d) 

^ = gi5C(3)A^/9g = 0.064736. (5.9e) 

Results in the relativistic cross-over region T ~ mth are plotted in Figs. ^ and ^. If one 
ignores the need to sum all ladder diagrams and only includes the one-loop diagram of 
Fig. I (after resumming self-energy corrections) then one underestimates the shear viscos- 
ity by roughly a factor of four. The analogous error for the bulk viscosity leads to an 
0(mpjjyg/A^T) result which scales completely incorrectly with A. 



^^The result ( |5.9| c) was not included in j^. In addition, the evaluation of T/^n in ref. 0| contained a 
numerical error which affected the plot of C shown in that paper. Recomputed values have been used in 



our Fig. I and Eq. (pjjp) 
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FIG. 4. Numerical results for the shear viscosity. The straight line shows the 0{T^/}?) asymptotic 
behavior of ?/. 
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FIG. 5. Numerical results for the bulk viscosity. The solid line shows the (mth/T)^ In^ (nrnth/T) 
behavior of Eq. (Bj 
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VI. CONCLUSIONS 



The analysis of this simple scalar field theory illustrates a number of points which are 
applicable to any relativistic field theory: 

a) The diagrammatic evaluation of transport coefficients is a remarkably inefficient ap- 
proach. An infinite set of rather complicated diagrams must be summed, merely to 
obtain the leading weak coupling behavior. 

b) The bulk viscosity depends on particle number changing processes and is sensitive 
to soft momenta, whereas the shear viscosity is determined by two body elastic 
scattering cross sections at typical momenta. The ratio of the bulk to the shear 
viscosities varies from very small (0(A^)) to exponentially large depending on the 
temperature. Hence, crude estimates such as C ~ ''7 which have appeared in 



the literature p,10| cannot generally be trusted. 



c) At high temperature, the existence of an effective kinetic theory adequate for com- 
puting transport coefficients depends crucially on the theory being weakly coupled, so 
that mean free paths are large compared to the wavelengths of relevant excitations. 
In, for example, high temperature QCD, it is unclear if the bulk viscosity can be 
correctly computed with any kinetic theory since the effective coupling of excitations 
with soft 0{g^T) momenta is not small. 

It is tempting to view the derivation of kinetic theory from the underlying field the- 
ory, and the derivation of the hydrodynamic constitutive equation (1.2) from the effective 
kinetic theory, as two different stages of a "real time renormalization group". At each 
stage, one is eliminating irrelevant degrees of freedom from the description of dynamics at 
successively lower frequency or momentum scales. We have little doubt that this notion of 
a real time renormalization group is essentially correct. However, we are unaware of any 
useful framework for defining a real time renormalization group which will systematically 
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transform the basic dynamical formulation from a quantum field theory to kinetic the- 
ory, or ultimately to classical hydrodynamics. In contrast to the situation for equilibrium 
Euclidean space observables [|11|], how to repackage the cumbersome diagrammatic analy- 
sis of in simple renormalization group terms is poorly understood. The diagrammatic 
treatment does not cleanly separate different frequency scales, as shown, for example, by 
the necessity of resumming both the real and imaginary parts of the on-shell single particle 
self energy in order to regulate individual cut diagrams, even though only the real part of 
the self energy appears explicitly in the resulting kinetic theory. The imaginary self energy, 
or single particle lifetime, should be viewed as an output of the effective kinetic theory, not 
an input parameter. A true real time renormalization group approach should allow one to 
derive completely the effective kinetic theory before treating any of the physics for which 
the kinetic theory description is adequate. Furthermore, a useful renormalization group 
framework should allow one to calculate corrections systematically, at least in weakly cou- 
pled theories. Although an effective kinetic theory did emerge in the analysis of the leading 
weak coupling behavior, it is unclear whether subleading corrections can be incorporated 
within a kinetic theory framework, since quantum coherence effects are only suppressed by 
a power of A. We hope that future investigations will shed light on some of these issues. 
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